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ABSTRACT 

We examine the accretion properties in a sample of 42 hard (3-60keV) X-ray 
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selected nearby broad-line AGNs. The energy range in the sample is harder than 
that usually used in the similar previous studies. These AGNs are mainly com- 
plied from the RXTE All Sky Survey (XSS), and complemented by the released 
INTEGRAL AGN catalog. The black hole masses, bolometric luminosities of 
\ AGN, and Eddington ratios are derived from their optical spectra in terms of the 

broad H/5 emission line. The tight correlation between the hard X-ray (3-20keV) 
and bolo metric/line luminosity is well identified in our sample. Also identified 
is a strong inverse Baldwin relationship of the H/3 emission line. In addition, all 

OO 

these hard X-ray AGNs are biased toward luminous objects with high Edding- 
ton ratio (mostly between 0.01 to 0.1) and low column density (< 10 22 cm -2 ), 
which is most likely due to the selection effect of the surveys. The hard X-ray 
luminosity is consequently found to be strongly correlated with the black hole 
mass. We believe the sample completeness will be improved in the next few years 
by the ongoing Swift and INTEGRAL missions, and by the next advanced mis- 
sions, such as NuSTAR, Simbol-X, and NeXT. Finally, the correlation between 
RFe (=optical Fe II/H/3) and disk temperature as assessed by T oc {L/ L-^^M^ 
leads us to suggest that the strength of the Fe II emission is mainly determined 
by the shape of the ionizing spectrum. 



Subject headings: galaxies: Seyferts — X-rays: galaxies — quasars: emission 
lines 
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1. Introduction 

It is now generally believed that the power of Active Galactic Nuclei (AGNs) is extracted 
through the accretion of gas onto a central supermassive black hole (SMBH). Such energy 
mechanism means that AGNs are characterized by their strong hard X-ray emission [hv > 
2keV), which is widely used as direct evidence suggesting the existence of a nuclear accretion 
activity (L 2 -8kcV > 10 42 er S s s_1 > e -S-> Silverman et al. 2005; Brusa et al. 2007; Hasinger 
et al. 2005). The commonly accepted model is that the hard X-ray emission from AGNs is 
primarily produced by the inverse Compton scattering of the UV/soft X-ray photons emitted 
from the accretion disk (e.g, Zdziarski et al. 1995, 2000; Haardt & Maraschi 1991; Kawaguchi 
et al. 2001). The absorption-corrected X-ray photon spectra within the energy band 2-10keV 
could be best described as a cut-off powerlaw with index T ~ 1.9 (e.g., Zdziarski et al. 1995; 
Reeves & Turner 2000; Piconcelli et al. 2005; Dadina 2008; Panessa et al. 2008). The 
synthesis spectra become flat beyond lOkeV because of the Compton reflection caused by 
the ionized surface of the accretion disk (e.g., George & Fabian 1991). 

Because the hard X-ray emission from central engine can penetrate the obscuration 
material much more easily than lower energy emission, it possesses particular importance in 
testing the traditional unified model (Antonucci 1993) in Seyfert 2 galaxies (e.g., Moran et al. 
2002; Cardamone et al. 2007). The studies of the Chandra and XMM-Newton observatories 
showed that the cosmic X-ray background (CXRB) at 2-30keV might be contributed by 
many unknown obscured AGNs which are predicted by the CXRB models (e.g., Comastri 
et al. 1995; De Luca & Molendi 2004; Worsley et al. 2005; Gilli et al. 2007; Severgnini et 
al. 2003; Levenson et al. 2006). Taking into account of the issue of co-evolution of AGN 
and bulge of its host galaxy (e.g., Heckman et al. 2004; Kauffmann et al. 2003; Wang et 
al. 2006; Wang & Wei 2008 and references therein), the hard X-ray emission from AGNs is 
also an important tool in detecting and separating AGN's contribution from circumnuclear 
star formation activity. The X-ray luminosities of known most X-ray luminous starforming 
and elliptical galaxies are not higher than L x = 10 42 erg s _1 (Zezas et al. 2003; Lira et al. 
2002a, 2002b; O'Sullivan et al. 2001). 

Heckman et al. (2005) identified a very tight correlation between the hard X-ray (3- 
20keV) and [O III] luminosities in a sample of hard X-ray selected AGNs when they performed 
a comparison between the hard X-ray selected and [O III] emission-line selected AGNs. In 2- 
lOkeV bandpass, similar correlations were identified in the Palomar optically selected AGNs 
by Panessa et al. (2006). On the contrary, a very weak L[om]-£x correlation was identified 
in the AGNs selected by their bright [O III] emission lines (Heckman et al. 2005). The result 
suggests that many AGNs might be missed in the hard X-ray survey. Moreover, Netzer et 
al. (2006) claimed that the L[om]/ix ratio depends on the X-ray luminosity. 
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The questions are therefore naturally raised: why is the L[om]-ix correlation broken in 
some kind of AGNs? which parameters (or what are the physical reasons that) determine the 
correlation? Are the hard X-ray selected AGNs particular in some parameters? Both black 
hole mass (Mrh) and Eddington ratio (L/L-^^) are two key parameters determining the 
observed properties of AGNs. In addition, with the development of the technology in hard 
X-ray detection, a major advance in studying AGN hard X-ray emission will be achieved in 
the next a few years due to the launch of new missions with enhanced hard X-ray capability 
(in sensitivity and imaging), such as NuSTAR, Simbol-X, and NeXT (e.g., Takahashi et al. 
2008; Ferrando et al. 2003). The study on the existent surveys certainly prepares the ground 
for the future surveys. 

In this paper, we examine the optical spectral properties of a sample of 42 hard X-ray 
selected broad-line AGNs, which allows us to investigate the properties of M BH and L/L Edd 
in these hard X-ray AGNs. The sample is mainly compiled from the AGN catalog of the 
RXTE 3-20keV All Sky Survey (XSS, Sazonov & Revnivtsev 2004), and complemented by 
the AGN catalog released by the INTEGRAL all-sky hard X-ray survey (Bassani et al. 2006a 
and references therein). Note that the energy range used in this work is extended to the 
Compton reflection region, and is harder than the range (i.e., 0.5-10 keV) usually used in the 
similar previous studies. The paper is organized as follows. Section 2 and 3 shows the sample 
selection and spectroscopic observations, respectively. The data reduction is described in the 
next section. Section 5 presents the results and discussions. The cosmology with h = 0.7, 
Vl M = 0.3, and f2 A = 0.7 (Bennett et al. 2003) is adopted in our calculations throughout the 
paper. 



2. The Sample 

Our sample is mainly compiled from the RXTE 3-20keV All-Sky Survey (XSS, Revnivt- 
sev et al. 2004). Sazonov & Revnivtsev (2004, hereafter SR04) identified 95 nearby AGNs 
with z mcdiim ~ 0.035 detected in XSS. The survey covers about 90% of the sky at |6| > 
10 ° . The sensitivity of the survey in the energy band 3-20 keV is better than 2.5 x 
10~ n ergs s _1 cm~ 2 . SR04 therefore provides a nearly complete sample for bright hard X-ray 
selected AGNs, although the sample is still biased against sources with absorption column 
N n > 10 23 cm- 2 . The IBIS telescope (Ubertini et al. 2003) onboard the INTEGRAL ob- 
servatory (Winkler et al. 2003) provides a better capability for heavily obscured objects 
because of its good sensitivity beyond 20 keV. After about three years INTEGRAL obser- 
vation, their all-sky survey project allows Bassani et al. (2006a) to identify 62 AGNs in the 
energy bandpass 20-100keV above a flux of 1.5 x 10~ n ergs s _1 cm~ 2 by analyzing 11,300 
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INTEGRAL points (see also in Bassani et al. 2006b; Bird et al. 2006; Beckmann et al. 
2006). 

Both catalogs are then combined to enlarge our sample content, and to attempt to 
alleviate the bias against the obscuration. We further limit the objects to the broad-line 
AGNs according to the identification provided by previous optical spectroscopy. The BAT 
instrument (Barthelmy et al. 2005) onboard the Swift satellite (Gehrels et al. 2004) provides 
an all-sky hard X-ray survey with a sensitivity down to a few xl0~ n ergs s _1 cm -2 in higher 
energy band 14-195 keV. Markwardt et al. (2005) identified 44 previous known AGNs in the 
Swift survey. However, considering the spectroscopic observation condition (see below), only 
three objects are not listed in our sample. It is worthy noting that the ongoing Swift/BAT 
and INTEGRAL survey could provide a larger sample of ~ 100 hard X-ray selected AGNs 
at the end of our program (i.e., Krivonos et al. 2006; Sazonov et al. 2007; Tueller et al. 
2008). We will perform a subsequent spectroscopic study on these AGNs in the next work. 

Figure 1 shows the column density (Nh) distribution for our sample. Nh is complied 
from SR04, and from Sazonov et al. (2007) for the objects detected by INTEGRAL only. 
Briefly, these authors determined Nh by fitting the available spectra from different X-ray 
instruments by an absorbed powerlaw. Values of Nh less than 10 22 cm -2 were ignored in the 
fitting of SR04 and Sazonov et al. (2007). As shown in the figure, the current sample is still 
strongly biased against AGNs with large column density. 



3. Spectroscopic Observations 

Both because of the constraint of the observatory site and because of the instrumental 
capability, the spectroscopic observations are only carried out for the bright (m v < 16.5mag) 
objects located in the northern sky with declination 5 > —20 ° . In total, there are 53 
objects fulfilling the selection criterion. Six of these objects are not observed because of 
the poor weather conditions. Among the remaining objects, five objects (i.e., XSSJ 02151- 
0033, XSS J11570+5514, XSS J22363-1230, NGC788, IGR J21247+5058) are excluded in 
subsequent spectral modeling because of the weakness or absence of broad Ha and/or H/3 
emission lineal 



1 Our spectroscopic observations show that three objects, i.e., XSS J02151-0033, XSS J11570+5514, and 
XSS J22363-1230, could be further classified as Seyfert 1.9-like galaxies without broad components of Balmer 
emission lines except a broad Ha emission. In addition, the continuum is dominated by the contribution 
from starlight of host galaxy rather than AGN in XSS J02151-0033 and XSS J11570+5514, which precludes 
the measurements of the Fe II complex. We include NGC 788 in our observation program because multi- 
classification is shown in the NED (see also in Beckmann et al. 2006). Our spectrum indicates a Seyfert 2 
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The total 42 high (or intermediate high) quality optical spectra were obtained by us- 
ing the National Astronomical Observatories, Chinese Academy of Science (NAOC), 2.16m 
telescope in Xinglong observatory during several observing runs carried out from November 
2005 to Mar 2008. The spectra were taken by the OMR spectrograph equipped with a back- 
illuminated SPEC10 1340x400 CCD as detector. A grating of 300g mm" 1 and a slit of 2.0" 
oriented in south-north direction were used in our observations. This setup results in a final 
spectral resolution ~ 9A as measured from both comparison spectra and night sky emission 
lines. The blazed wavelength was fixed to be 6000A, which provides a wavelength coverage 
of 3800-8300A in observer frame. This attempt covers both Ha and Hf3 region in most 
spectra because of their small redshifts. Each object was observed successively twice. The 
two exposures were combined prior to extraction to enhance the S/N ratio and eliminate the 
contamination of cosmic-rays easily. The exposure time for each frame is generally between 
900 and 3600 seconds depending on the brightness of the object. The wavelength calibration 
associated with each object was carried out by the helium-neon-argon comparison arcs taken 
between the two successive frames. The arcs were obtained at the position being nearly iden- 
tical to that of the specified object. Two or three Kitt Peak National Observatory (KPNO) 
standard stars (Massey et al. 1988) were observed per night for both flux calibration and 
removal of the atmospheric absorption features. All the objects were observed as close to 
meridian as possible. Table 1 shows the log of observations of the 42 objects listed in our 
sample. 



4. Data Reduction and Emission Line Measurements 

The standard procedures using the IRAF package @ are adopted by us to reduce the 
unprocessed frames. The CCD reductions include bias subtraction, flat-field correction, 
and cosmic-rays removal before the extraction of signal. One- dimensional sky-subtracted 
spectra are then wavelength and flux calibrated. The uncertainties of the wavelength and 
flux calibrations are no more than lA and 20%, respectively. The Galactic extinctions 



nucleus in the object. In fact, the object is classified as a Seyfert 2 galaxy with polarized broad emission 
lines in the catalogue of quasars and active nuclei (Veron-Cetty & Veron 2006). Our spectrum taken in Dec 
20, 2006 shows that the spectrum of IGR J21247+5058 is dominated by a typical late type star, which is 
consistent with Masetti et al. (2004). By identifying the broad emission bump around 6700A as Ha, these 
author argued that the background AGN is aligned with a F- or early G-type star by chance. 

2 IRAF is distributed by the National Optical Astronomical Observatories, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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are corrected by the color excess, the parameter E(B — V) taken from the NASA/IAPC 
Extragalactic Database (NED), assuming an Ry = 3.1 extinction law of Milk Way (Cardelli 
et al. 1989). The spectra are then transformed to the rest frame, along with a /^-correction, 
according to the narrow peak of 

4.1. Fe II subtractions 

In order to determine the strength of the optical Fe II complex and to reliably profile the 
emission lines in H(3 region, the Fe II complex should be modeled and then removed from each 
spectrum at first. We simply adopt the empirical template of the Fe II complex introduced 
in Boroson & Green (1992, hereafter BG92). For each object, the template is first broadened 
to the FWHM of the H(3 broad component by convolving with a Gaussian profile (BG92). 
The broadened template and a broken powerlaw are combined to model the continuum by 
a x 2 minimization for individual object. The minimization is performed over the rest frame 
wavelength range from 4300 to 7000A, except for the regions around the strong emission lines 
(i.e., Ha+[N II], [S II], R(3, FFy, [O III], He II). The fittings are illustrated in the left column 
of Figure 2 for three typical cases. The flux of the Fe II complex is integrated between the 
rest wavelength 4434 and 4686A. The Fe II complex is not measured in XSS J11067+7234 
because its continuum is dominated by the starlight from the host galaxy. 

4.2. Emission lines measurements 

After removing the Fe II blends and continuum, the isolated AGN emission lines at H/3 
region are modeled by the SPECFIT task (Kriss 1994) in IRAF package. In each object, 
each of the [O III] AA4959, 5007 doublet is modeled by a single Gaussian component. The 
flux ratio of the doublet is fixed to be 3 (e.g., Dimitrijevic et al. 2007). In principle, each 
H(3 profile is modeled by a set of two Gaussian components: a narrow and a broad Gaussian 
component (denoted as H/3 n and H/5b, respectively). The profile modelings are illustrated 
for the three cases in the right column in Figure 2. The deblendings of the H/5 emission line 
are, however, difficult in 15 objects because no obvious transition between the narrow and 
broad components is observed. In these 15 cases, the reported H/3 b measurements include 
both components, and the contribution of Hj3 Q is usually expected to be less than 3% for 
typical broad-line AGNs as suggested by BG92. 

The described modeling works very well in most objects, except in three ones. In 
XSSJ22539-1735, a very broad H/5 component (FWHM ~ 10,000 km s" 1 ) is required to 
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properly reproduce the observed profile. Such component has been reported in H/5 and 
highly ionized emission lines in previous studies (e.g., Sulentic et al. 2000b; Wang et al. 2005; 
Veron-Cetty et al. 2007; Marziani et al. 2003; Mullaney & Ward 2008; Hu et al. 2008). A 
double peaked and a boxy H/3 profile is observed in XSS J18408+7947 and XSS J23073+0447, 
respectively. In all of the three cases, each H/3 line is modeled by three Gaussian components, 
i.e., a narrow and two broad components. The modeled narrow component is then subtracted 
from the observed spectrum to derive a residual profile which is then used to measure FWHM 
and integrated flux. 



5. Results and Discussion 

Table 2 lists the following items measured from the spectra. Column (2) and (3) lists the 
equivalent widths (EWs) of H/3 n and H/3b, respectively. All the EWs refer to the continuum 
level determined from the continuum fitting at the position of wavelength A4861. The line 
ratio RFe (=Fe II/H/^) is shown in column (4). Column (5) lists the luminosity of the 
H/3 broad component, and Column (6) the FWHM of the broad H/3 in units of km s _1 . 
All quoted widths are not corrected for the intrinsic resolution, since the widths are much 
broader than the resolution. Column (7) shows the luminosity of the [O III] emission line. 
Both [O III] and H/5b luminosities are calculated from the cosmology-corrected luminosity 
distance [f| given by NED. The next two columns show the absorption-corrected hard X-ray 
luminosities in the bandpass 3-20 keV (XSS) and in 17-60 keV (INTEGRAL). 

Six out of the 42 AGNs are detected by INTEGRAL only. In order to estimate their 
X-ray luminosities in the bandpass 3-20 keV from the luminosities in 17-60 keV, a trans- 
formation is derived in terms of the 18 objects detected by both surveys. The left panel of 
Figure 3 shows the relationship between the two sets of luminosity. A unweighted fitting 
yields a relationship logLi7_ 60 k c v = 1-03 log L 3 _ 2 okcV — 1-39 with a standard deviation of 0.25. 
After transform the luminosities in the bandpass 17-60 keV to 3-20 keV for the six objects, 
the strong correlation between the [O III] and 3-20 keV luminosity is shown in the right panel 
of Figure 3 for the total 42 objects. The correlation is highly consistent with that derived 
in Heckman et al. (2005), which firmly demonstrates the accuracy of our observations and 
calibrations. 



3 Because the objects listed in our sample are generally very nearby, the radial velocity is corrected by a 
redshift of 0.017877 to the reference frame defined by the 3K Microwave Background Radiation. 
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5.1. Hard X-ray vs. bolometric luminosity 

The detection of hard X-ray emission from a nucleus is regarded as strong evidence of 
accretion activity occurring around the central SMBH. The correlations between X-ray lumi- 
nosity and luminosities of optical emission lines (e.g., Ha, H/3, [O III]) have been extensively 
established in the previous studies (e.g., Elvis et al. 1984; Ward et al. 1988; Mulchaey et 
al. 1994; Heckman et al. 2005; Panessa et al. 2006). However, the correlations seem to 
depend on sample selection (Heckman et al. 2005) and on X-ray luminosity itself (Netzer et 
al. 2006). In this paper, we estimate the luminosity of AGN at rest wavelength 5100A from 
the H/?;, component according to the calibration given by Greene & Ho (2005): 

/ Ijw \ °- 883 

L 5100 = 7.31xl0 43 ( ^ ) ergss- 1 (1) 
M0 42 erg s i / 

The bolometric luminosity is then derived by adopting the usually used bolometric correction 
L ho \ Rj 9AL A (5100A) (e.g., Kaspi et al. 2000). The significant correlation between L bol and 
-^3-20keV is shown in the left panel in Figure 4. A unweighted least square fitting gives a 
relationship: 

log L x = (0.91 ± 0.06) log L bo i + (3.04 ± 2.78) (2) 

or Lboi/Lx ~ 10 for simplification. The correlation strongly indicates a close link between 
the hard X-ray emission and ionizing radiation emitted from the accretion disk. Note that 
the bolometric luminosity spans four orders of magnitude, down to L ho \ ~ 10 43 ergs s _1 , 
which is close to the formal definition of low luminous AGNs (i.e., L DO i < 10 43 ergs s _1 , Ho 
2003). These faint galaxies generally radiate at low state with Eddington ratio marginally 
exceeds 0.01. 



5.2. Black Hole Mass and Eddington Ratio 

It is now generally believed that black hole mass (Mbh) and specific accretion rate 
(L/L-Edd) are two basic parameters determining the properties of AGNs. The role of L/L-^dd 
in driving the Eigenvector I (EI) space has been extensively investigated in numerous 
previous studies (e.g., Boroson 2002; Sulentic et al. 2006; Xu et al. 2003; Grupe 2004; 



4 The EI space is one of the key properties of AGN phenomena. It was first established by BG92 who 
analyzed the optical spectra of 87 bright PG quasars. In addition to the anti-correlation between the intensity 
of Fe II and [O III], the EI space has been subsequently extended to ultraviolet and soft X-ray bands (i.e., 
FWHM(C IV) and T s , e.g., Wang et al. 1996; Xu et al. 2003; Grupe 2004; Sulentic et al. 2007; see Sulentic 
et al. 2000a for a review). 
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Zamanov & Marziani 2002; Marziani et al. 2001), because of the great progress in the 
calibration of the -Rblr — L relationship (e.g, Kaspi et al. 2000, 2005, 2007; McLure & 
Jarvis 2004; Vestergaard & Peterson 2006; Peterson et al. 2004; Bentz et al. 2006) due to 
the recent great advance in the reverberation mapping (e.g., Kapsi et al. 2000; Peterson & 
Bentz 2006). We refer the reads to McGill et al. (2008, and references therein) for a summary 
of the existing formula used to calculate Mbh basing upon "single-epoch" observation. In 
this work, the black hole mass is estimated by the width and luminosity of Hfo in individual 
object, according to the scaling law obtained by Greene & Ho (2005). 



For each object, the estimated Mbh and L/L Edd is listed in Column (10) and (11) in Table 
2, respectively. The bolometric luminosity L bol is estimated from the H/3b component as 
described above. 

The calculations of both Lboi and L/L Edd allow us to find that the current sample 
is strongly biased against sub-luminous AGNs usually with low L/L Edd . In fact, the left 
panel of Figure 4 shows the lack of AGNs with Lboi < 10 43 ergs s -1 . Figure 5 displays the 
distributions of the L/L Edd (left panel) and Mbh (middle panel) for the total 42 AGNs. As 
shown in the middle panel, the Mbh is sampled within a wide range (~ 3 dex) from 10 6 M Q 
to 1O 9 M with a peak at 1O 7 ~ 8 M . In contrast, the left panel shows that L/L Edd distributes 
in a relatively narrow range as compared with the previous studies. The total range of our 
L/L Edd spans from 0.01 to 1. In particular, about ~ 60% objects listed in the sample have 
L/L Edd between 0.01 and 0.1. However, the L/L Edd of bright local AGNs usually spans at 
least three orders of magnitude from 1 to 0.001 (e.g., Woo & Urry 2002; Boroson 2002). 
Basing upon the 2-10keV X-ray luminosity, Panessa et al. (2006) indicated that the L/L Edd 
of the Palomar optically selected AGNs ranges from 10~ 5 to 0.1. As an additional test, the 
right panel of Figure 4 shows the M B h vs. L 3 _ 2 okcV plot. L 3 _ 2 okcV as a function of M B h is 
over-plotted as dashed lines for three different L/L Edd (i-e., 1, 0.1, and 0.01). As shown in 
the plot, a majority of our objects are located below the line with L/L E dd = 0.1 and above 
the line with L/L Edd = 0.01. 

In summary, the hard X-ray selected AGNs listed in our sample are luminous AGNs 
with a wide range of Mbh but a nearly constant accretion activity (i.e., L/L Edd ), which is 
likely due to the selection effect of the survey that is biased toward X-ray luminous objects. 
In active (or luminous) AGN, the optical/UV ionizing radiation is believed to be emitted 
from a standard geometric thin disk (e.g., Shakura & Sunyaev 1973). The low energy photon 
is comptonized by a hot corona to produce the hard X-ray emission below lOkeV (Haardt 
& Maraschi 1991). A Compton recoil of the soft photon is required to occur on the ionized 




(3) 
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surface of the accretion disk to produce the emission spectrum beyond lOkeV. The main 
observation feature of the reflection is a bump peaked at about 30keV (e.g., George & 
Fabian 1991; Zycki et al. 1994; Ross & Fabian 2002, 2005). 

The bias towards active AGNs could be possibly caused by the fact that either intensive 
Compton reflection takes place only in AGNs at high state with large L/L E dd or the surveys 
are biased against the Compton-thick objects. In the first case, the theory of the Compton 
reflection predicts that the X-ray emission contributed by the reflection depends primarily 
on the X-ray ionizing parameter, and secondarily on the UV radiation produced by the 
dissipation inside the accretion disk. In the second case, our analysis implies a possible 
connection between the less X-ray absorption and high L/L-^dd i n broad-line AGNs. In fact, 
we selected the objects with regardless of their X-ray spectral properties. Figure 1 shows 
the lack of objects with large column density in the sample. We believe that the sample 
completeness would be improved by including the ongoing Swift/BAT survey with large 
effective collecting area and harder energy bandpass (14-195 keV) in the future studies. 

Figure 4 shows that the hard X-ray luminosity is strongly correlated with the M BH 
in our sample. In fact, the correlation is naturally expected given the tight L3_20k e v vs- 
Lboi correlation and nearly constant L/LEdd- The Lx vs. Mbh correlation provides us a 
potential estimate of black hole mass for luminous AGNs within 0.01 < L/LEdd < 0.1. The 
following relationship is obtained by us through a least square fitting: \og(M-Bii/M & ) = 
(0.82 ± 0.09) logL 3 _20kcV - (28.23 ± 3.78). Our results conflicts with Panessa et al. (2006) 
and Pellegrini (2005) who did not find the correlation, but is in agreement with Kiuchi et al. 
(2006). Panessa et al. (2006) investigated a sample of 47 nearby Seyfert galaxies selected 
from the Palomar spectroscopy (Ho et al. 1997). The luminosity obtained by different X-ray 
instruments is down to L 2 _iokcV ~ 10 38 ergs s _1 . The studies in Pellegrini (2005) are based 
on the Chandra observations down to a luminosity L 2 _iokcV ~ 10 38 ergs s _1 . Kiuchi et al. 
(2006) used the broad-line AGN sample detected by the ASCA Large Sky Survey (ALSS) 
and ASCA Medium Sensitivity Survey in the northern sky (AMSSn) with a detection limit 
of a few xl0~ 13 ergs s _1 cm -2 in 2-10keV bandpass. As discussed above, it is worthy noting 
that the correlation most likely does not reflect the physics but a selection effect of the 
surveys. 



5.3. Fe II ratio vs. disk temperature relation 

As a key parameter in the EI space, the Fe II ratio (RFe) is defined as the ratio of the 
optical Fe II complex to Hf3. Although the total Fe II emission increases by fourfold with 
reasonable ionization parameter in AGN (Korista et al. 1997), traditional photoionization 
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models can not explain the strong Fe II emission in optical and UV bands (e.g, Netzer & 
Wills 1983; Joly 1987; Collin-Souffrin et al. 1988). At present, the problem is only slightly 
alleviated by the major improvements in the atomic data and by the improved treatment 
of the line excitation process (Sigut & Pradhan 2003; Baldwin et al. 2004). We refer 
the reads to Collin & Joly (2000) for a summary of the mechanisms that can enhance the 
Fe II emission. BG92 put forward a picture in which the RFe is determined by the vertical 
structure of the accretion disk. The vertical structure is governed by L/L Edd . A large L/L Edd 
leads to a large X-ray heated volume that could generate large Fe II emission. Sulentic et 
al. (2000b) developed a semi-analysis model in which the RFe depends on L/ L Edd as RFe oc 

0. 5510g(L/L Edd ). Although the RFe is found to generally increase with L/L Edd in optical 
bright quasars, Netzer et al. (2004) found that a number of high-z quasars deviate the trend, 

1. e., with very small RFe but large L/L Edd (see also in Netzer & Trakhtenbrot 2007). Using 
the large database provided by the Sloan Digital Sky Survey, Netzer & Trakhtenbrot (2007) 
recently suggested that the enhanced RFe is mainly caused by increased metal abundance. 

In the current sample, the distribution of RFe is shown in the right panel of Figure 5. RFe 
uniformly ranges from 10~ 3 to 1. On the contrary, L/L Edd distributes in a quite narrow range 
as describe above. In fact, no correlation between RFe and L/L Edd is found in our sample 
(see the upper panel in Figure 6), which motivates us to suspect that RFe does not depend on 
L/L Edd only. RFe is plotted against the characteristic disk temperature T max in the bottom 
panel in Figure 6. The temperature scales with L/L Edd and Mrh as predicted by the standard 
geometric thin disk model (e.g., Shakura & Sunyaev 1973). The exact formula of disk 
temperature depends on various accretion disk models. For a rapidly rotating Kerr hole, with 
a spin parameter a* = 0.998 and efficiency of 0.31, we have T max = 10 5 56 (L/L Edd ) 1 / 4 M BI ! [ ^ 4 K. 
The diagram indicates an obvious correlation between the two parameters. A spearman rank- 
order test calculated by survival analysis yields a formal correlation coefficient r s = 0.414 
(P = 0.0088, where P is the probability of null correlation). The correlation is not highly 
significant probably because of the small sample size. The estimated temperature spans 
about one order of magnitude (AlogT max ps 1-25), corresponding to a factor of ~ 18. The 
correlation then suggests a trend of more intensive Fe II emission for higher disk temperature. 
A marginal dependence of the continuum shape of QSOs on T max was recently identified by 
Bonning et al. (2007) who compared the observations of SDSS with the NLTE models of 
accretion disk. The current result means that the strength of the Fe II emission is likely 
controlled by the spectral shape of the ionizing continuum. 
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5.4. Equivalent width of H/3 vs. Lx 

The physical reason of the absence of the Baldwin relationship (Baldwin 1977) for low 
ionization emission lines is still an open question. In fact, a weak inverse Baldwin relationship 
for H/3 has been demonstrated by recent studies basing upon large AGN samples (e.g., Croom 
et al. 2002; Greene & Ho 2005). We identify a tight, positive correlation between L x and 
EW(H/3 b ) in our hard X-ray selected AGNs (i.e., an inverse Baldwin relationship). Figure 
7 presents the correlation with correlation coefficient r = 0.611 (P = 10~ 4 ) estimated by 
the spearman rank-order analysis. EW(H/3&) roughly scales with hard X-ray luminosity as 
EW(H/3 b ) oc L^ 39 . Wilkes et al. (1999) identified a marginal Baldwin effect in H/3 line. 
Noted that they examined only the luminous local quasars with Li-iokcV > 10 44 ergs s -1 . 

Although many models are developed to explain the Baldwin effect for high ionization 
emission lines (e.g., C IV, Wandel 1999; Korista et al. 1998; Shields et al. 1995; Wills 
et al. 1999; Baskin & Loar 2004; Bachev et al. 2004), these models can not explain the 
difference between C IV and H/3. Croom et al. (2002) suggested that the inverse Baldwin 
relationship could be explained if the longer wavelength continuum contains emission from 
other components (e.g., thermal dust emission, non-thermal radio emission, starlight). We 
estimate the possible contribution of the unknown sources as follows. We start from the 
relationship logL x = 39.08 + 2.56 log EW(H/3), and re-write EW(H/3) = L{H(3)/{L' + L c ), 
where L c and V is the AGN luminosity and luminosity of other unknown sources at the H/3 
wavelength, respectively. The X-ray luminosity could be replaced by L c given Eq. (2) 
and the bolometric correction factor of 9. Replacing L(Hj3) as L c given Eq.(l) finally 
yields a relationship log(l + L'/L c ) = 9.48 — 0.211ogL c . Considering the typical case with 
L c ~ 10 44 erg s -1 , about 40% of the observed continuum at H/3 wavelength is estimated to 
be contributed by the unknown sources. 



6. Conclusion 

The properties (L/LEdd and Mbh) of accretion onto SMBH are examined in a sample of 
42 hard X-ray selected (3-60keV) broad-line AGNs in terms of their optical spectra taken by 
us. The energy range is harder than that usually used in the similar previous studies. These 
AGNs are mainly compiled from the RXTE All Sky Survey (Sazonov & Revnivtsev 2004), 
and are complemented by the released INTEGRAL AGN sample (Bassani et al. 2006a). 
The statistical analysis allows us to draw the following conclusions: 

1. We confirm the tight correlation between the hard X-ray and optical emission line 
luminosities (and bolometric luminosity) in our sample, which suggests a close link 
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between the hard X-ray emission reflected by the ionized surface of the accretion disk 
and UV/optical radiation. Using the hard X-ray luminosity, a strong inverse Baldwin 
relationship of the H/3 emission line is identified in the sample. 

2. The hard X-ray selected broad-line AGNs listed in the sample are found to be strongly 
biased toward luminous AGNs with high Lj L Edd and low column density. Since Lj L Edd 
is constant (mostly between 0.01 and 0.1) in a first order approximation, the hard X- 
ray luminosity is strongly correlated with the black hole mass in our sample, which is 
most likely due to the selection effect of the surveys. 

3. Although the RFe parameter is independent on L/L Edd in our sample, it is found to 
be correlated with the accretion disk temperature as assessed by T oc (L/L E dd)^BH- 
This result implies that the strength of the Fe II emission is determined by the shape 
of the ionizing spectrum. 

Finally, it should be mentioned that a new era in AGN hard X-ray study will be opened 
in next a few years due to the launch of new missions with enhanced hard X-ray detection 
capability in not only sensitivity, but also imaging, such as Simbol-X, NeXT and NuSTAR. 
These missions will provide larger, and more complete samples to study the present open 
issues. 

We would like to thank the anonymous referee for his/her valuable comments that 
help to improve the paper. The authors are grateful to Todd A. Boroson and Richard 
F. Green for providing us the Fe II template. Special thanks go to the staff at Xinglong 
observatory as a part of National Astronomical Observatories, China Academy of Science 
for their instrumental and observational help. This search is supported by the NFS of China 
under grant 10503005. 
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Fig. 1. — Observed X-ray absorption distribution of the 42 hard X-ray selected broad-line 
AGNs listed in the sample. 
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Fig. 2. — left column: the modeling of the continuum and Fe II complex around B./3 region 
for three typical cases. In each panel, the blue curve shows the observed spectrum, and the 
overlaid red one the modeled spectrum (continuum + Fe II complex). The strong emission 
lines, especially the two strong Fe II complex around the H/3, are marked on the figure, right 
column: the modeling of the emission-line profile by a set of Gaussian components for the 
H/9 and [O III] doublet in the three cases, after subtract the modeled continuum and Fe II 
complex. In each panel, the observed profile is shown by the black curve. The three red 
short dashed lines show the modeled H/3 broad component and [O III] doublet. The modeled 
narrow H/? emission is presented by the long dashed line. 
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Fig. 3. — left panel: the luminosity in the energy bandpass 3-20 keV plotted against the 
luminosity in 17-60 keV for the 18 objects detected by both XSS and INTEGRAL. The 
solid line shows the best fit to the data without taking into account of the errors. The la 
dispersion is marked by the two dashed lines, right panel: the tight correlation between the 
[O III] emission line and 3-20 keV luminosity. 
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Fig. 4. — left panel: L 3 _ 2 okcV plotted against the bolometric luminosity estimated from the 
H/3 broad component. The solid line shows the best fit to the data: logL^ = (0.91 ± 
0.06) logLboi + (3.04 ±2.78). right panel: Distribution on the plot of Lx vs. black hole mass. 
The three dashed lines show the Lx as a function of the black hole mass for different L / LecM 
=0.01, 0.1, and 1. 
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Fig. 5. — Distribution of the Eddington ratio (left panel), black hole mass (middle panel), 
and parameter RFe (right panel). 
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Fig. 6. — top panel: RFe plotted against L/L Edd (r s = 0.237, P = 0.134). bottom panel: 
the same as the top panel, but for T max = 10 5 - 56 (L/L Edd ) 1 / 4 M BH ~ 1/4 - The diag ram shows a 
relatively strong correlation between the two parameters (r s = 0.414, P = 0.0088). In both 
panels, the points with zero flux of the Fe II emission are indicated by superposed arrows. 
The correlation coefficients are calculated by survival analysis. 
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Fig. 7. — Hard X-ray luminosity is plotted as a function of equivalent width of the H/3 
broad component. A spearman rank-order test yields a correlation coefficient r = 0.611 
(P = 10" 4 ). 
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Table 1: Log of spectroscopic observation 



Name 


Other name 




R.A. 


Dec] 


.. 


z 


Date 
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(s) 


m v 


AGN Type 


(1) 


(2) 




(3) 




(4) 




(5) 
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00 


36 


21.0 


45 


39 


54 


0.048 
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15 


15.45 
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02 


44 
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62 


28 


07 
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16 
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56 
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NLS1 


XSS J13530+6916 


MRK 279 


13 


53 


03.4 


69 


18 


30 


0.031 


1200 


2007 


Feb 


16 


14.57 


SI 


XSS J14181+2514 


NGC 5548 


14 


17 


59.5 


25 


08 


12 


0.017 


1800 


2008 


Mar 


03 


13.3 


SI 


XSS J15042+1046 


MRK 841 


15 


04 


01.2 


10 


26 


16 


0.036 


1200 


2007 


Feb 


16 


14.0 


SI 


XSS J15348+5750 


MRK 290 


15 


35 


52.3 


57 


54 


09 


0.030 


2400 


2008 


Mar 


13 


15.1 


SI 


XSS J17276-1359 


PDS456 


17 


28 


19.8 


-14 


15 


56 


0.184 


1800 


2007 


Sep 


07 


14.69 


RQQ 


XSS J17413+1851 


4C +18.51 


17 


42 


07.0 


18 


27 


21 


0.19 


2400 


2007 


Sep 


07 


16.43 


RLQ 


XSS J18196+6454 


H 1821+643 


18 


21 


57.3 


64 


20 


36 


0.297 


1800 


2006 


Aug 


16 


14.1 


RQQ 


XSS J18348+3238 


3C 382 


18 


35 


02.1 


32 


41 


50 


0.059 


1800 


2006 


Aug 


16 


15.5 


BLRG 


XSS J18408+7947 


3C 390.3 


18 


42 


09.0 


79 


46 


17 


0.056 


1800 


2006 


Aug 


15 


14.37 


BLRG 


1H 1934-063A 




19 


37 


33.0 


-06 


13 


05 


0.01059 


1800 


2006 


Aug 


16 


14.09 


SI 


NGC 6814 




19 


42 


40.6 


-10 


19 


25 


0.00521 


1800 


2007 


Nov 


04 


12.06 


SI. 5 


XSS J20404+7521 


4C +74.26 


20 


42 


37.3 


75 


08 


02 


0.1 


1800 


2006 


Sep 


15 


15.13 


RLQ 


XSS J20441-1042 


MRK 509 


20 


44 


09.7 


-10 


43 


25 


0.034 


1800 


2006 


Sep 


15 


13.0 


SI 


XSS J21128+8216 


S5 2116+81 


21 


14 


01.2 


82 


04 


48 


0.084 


2400 


2006 


Aug 


15 


15.7 


BLRG 


XSS J22423+2958 


AKN 564 


22 


42 


39.3 


29 


43 


31 


0.025 


1800 


2006 


Nov 


15 


14.55 


NLS1 


XSS J22539-1735 


MR2251-178 


22 


54 


05.8 


-17 


34 


55 


0.064 


1800 


2006 


Nov 


15 


14.36 


RQQ 


XSS J23033+0858 


NGC 7469 


23 


03 


15.6 


08 


52 


26 


0.016 


1800 


2006 


Aug 


15 


13.0 


SI 


XSS J23040-0834 


MRK 926 


23 


04 


43.5 


-08 


41 


09 


0.047 


1800 


2007 


Nov 


03 


14.6 


SI 


XSS J23073+0447 


PG 2304+042 


23 


07 


02.9 


04 


32 


57 


0.042 


2400 


2007 


Nov 


02 


15.44 


SI 



Table 2. List of properties of the hard X-ray selected AGNs 



Name EW(H/3 n ) EW(H/3 b ) RFc log[L(H/3 b )} FWHM(H/3 6 ) log[L([OIII])] log[L 2 _ 3 okev] log[L 17 _ 60k ev] log(M B H/M Q ) L/L Edd 

A A ergs s _1 km s -1 ergs s _1 ergs s _1 ergs s _1 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



XSSJ 00368+4557 




82.6 


0.35 


42.1 


2555 


41.3 


43.48 




7.42 


0.23 


QSO 0241+622 




43.6 


0.26 


42.1 


9401 


41.8 


44.34 


44.28 


8.54 


0.02 


3C111 




119.9 


0.00 


43.3 


3445 


42.8 


44.53 


44.59 


8.33 


0.30 


XSS J04331+0520 


10.5 


82.7 


0.24 


42.2 


5367 


42.1 


44.16 




8.11 


0.06 


XSS J05103+1640 


9.0 


74.5 


> 0. 


41.3 


4855 


41.5 


43.47 


43.58 


7.55 


0.04 


XSS J05162-0008 


2.3 


72.7 


0.55 


42.5 


6595 


41.5 


44. 


44.26 


8.48 


0.05 


XSS J05552+4617 


53.5 


103.2 


0.30 


41.2 


5761 


41.7 


43.56 


43.67 


7.60 


0.02 


Mark 6 


6.0 


37.4 


0.42 


41.4 


5627 


41.7 


43.32 


43.41 


7.74 


0.03 


XSS J07434+4945 


8.6 


65.8 


0.40 


41.6 


3815 


41.7 


43.45 




7.47 


0.07 


XSS J081 17+7600 




74.7 


0.52 


43.4 


3146 


42.3 


44.47 


44.43 


8.32 


0.39 


XSS J09204+1608 


6.9 


91.8 


0.29 


41.9 


6421 


41.6 


43.65 




8.11 


0.03 


XSS J09261+5204 




118.6 


0.08 


42.1 


1798 


41.8 


43.81 


44.17 


7.12 


0.47 1 


XSS J10231+1950 


5.1 


39.3 


0.007 


40.5 


4556 


40.6 


42.4 


42.67 


7.01 


o.oi^ 


XSS J11067+7234 




31.9 




41.1 


5020 


40.7 


43.04 




7.46 


0.03 


XSS J11417+5910 


7.3 


73.8 


0.13 


42.0 


4608 


41.9 


43.82 




7.91 


0.07 


XSS J12032+4424 




21.5 


0.61 


39.9 


1202 


40.0 


42.04 


41.87 


5.53 


0.20 


XSS J12106+3927 


37.2 


100.0 


0.12 


40.8 


6354 


41.3 


43.17 


43.37 


7.49 


0.01 


XSS J12164+1427 




120.3 


0.40 


42.5 


1804 


41.7 


44.22 




7.32 


0.60 


NGC4253 


11.3 


32.1 


0.72 


40.8 


2153 


41.3 


42.59 


42.7 


6.57 


0.13 


XSS J12206+7509 




95.6 


0.21 


42.5 


3037 


41.8 


43.88 




7.83 


0.23 


XSS J12288+0200 




71.3 


0.51 


44.1 


2970 


43.2 


45.91 


45.92 


8.66 


0.74 


XSS J12408-0516 




49.8 


0.54 


41.0 


3654 


40.4 


42.57 


43.04 


7.12 


0.05 


XSS J13420-1432 


15.3 


64.6 


0.66 


41.8 


4112 


41.8 


43.57 




7.67 


0.07 


XSS J13530+6916 


5.7 


59.6 


0.44 


42.2 


5147 


41.8 


43.94 




8.07 


0.06 


XSS J14181+2514 


7.8 


62.7 


0.04 


41.3 


10312 


41.5 


43.66 


43.08 


8.20 


0.008 


XSS J15042+1046 


9.4 


129.6 


0.33 


42.5 


5709 


42.2 


43.67 




8.33 


0.061 


XSS J15348+5750 


4.2 


23.5 


0.71 


41.1 


4740 


41.6 


43.43 




7.41 


0.03 


XSS J17276-1359 




74.6 


3.58 


44.1 


2957 


< 42.7 


44.93 




8.66 


0.75 


XSS J17413+1851 


1.8 


115.5 


0.04 


43.0 


13506 


42.2 


45.45 




9.35 


0.02 


XSS J18196+6454 


2.9 


65.4 


0.18 


44.3 


5797 


43.7 


45.61 




9.36 


0.23 


XSS J18348+3238 




111.4 


0.19 


42.8 


8340 


41.8 


44.58 


44.45 


8.84 


0.04 


XSS J18408+7947 


3.2 


118.8 


0.15 


42.5 


13799 


42.0 


44.39 


44.61 


9.13 


0.01 



Table 2 — Continued 



Name 


hj W (Xl/Jn ) 




lire 




r VVrilVl^rlpf) ) 


log[_H[Ulli]J] 


log[L2_30keVj 


log[Li7_60k c vJ 


lOg(iWBH/JM0 ) 


^/-^Edd 




X 
J V 


X 




ergs s ^ 


km s ^ 


ergs s ^ 


ergs s ^ 


ergs s ^ 










(3) 








(<) 


(°) 


1,9) 


(1UJ 




1H 1934-063A 




55.2 


0.94 


41.1 


1354 


41.0 


42 AS 


42.59 


6.29 


0.38 


NGC6814 




38.8 


0.75 


40.2 


5182 


40.2 


42.35 


42.47 


6.95 


0.01 


XSS J20404+7521 




85.4 


> 0. 


43.6 


9418 


42.5 


44.73 


44.98 


9.37 


0.05 


XSS J20441-1042 


10.7 


111.8 


0.21 


42.4 


2261 


42.1 


44.2 


44.13 


7.47 


0.36 


XSS J21128+8216 


6.3 


150.0 


0.14 


42.6 


4123 


42.3 


44.25 


44.65 


8.12 


0.13 


XSS J22423+2958 




56.6 


0.74 


41.7 


3649 


41.7 


43.42 




7.52 


0.09 


XSS J22539-1735 


5.7 


132.5 


0.02 


42.9 


4093 


42.6 


44.63 


44.62 


8.28 


0.16 


XSS J23033+0858 


13.7 


45.3 


0.30 


41.5 


2819 


41.6 


43.27 


43.40 


7.17 


0.12 


XSS J23040-0834 


3.4 


65.1 


> 0. 


42.5 


7083 


42.2 


44.23 


44.21 


8.51 


0.04 


XSS J23073+0447 


2.7 


42.0 


> 0. 


41.2 


7322 


41.2 


43.67 




7.85 


0.01 
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